Lower critical field Hd and barriers for vortex entry in Bi2Sr2CaCu208+<5 crystals 
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The penetration field Hp of Bi2Sr2CaCu2 08+i crystals is determined from magnetization curves 
for different field sweep rates dH/dt and temperatures. The obtained results are consistent with 
theoretical reports in the literature about vortex creep over surface and geometrical barriers. The 
frequently observed low-temperature upturn of Hp is shown to be related to metastable configura- 
tions due to barriers for vortex entry. Data of the true lower critical field Hci are presented. The 
low-temperature dependence of H^i is consistent with a superconducting state with nodes in the 
gap function. [PACS numbers: 74.25.Bt, 74.60.Ec, 74.60. Ge, 74.72.Hs] 



High-temperature superconductors (HTSC's) are 
strongly type II and as such their H — T phase dia- 
gram involves complete flux expulsion (Meissner state) 
at low fields and magnetic field penetration in the form 
of quantized flux-lines or vortices (mixed state) at higher 
fields. The onset of the mixed state, where vortex pene- 
tration becomes energetically favorable, is defined as the 
lower critical field Hd oc 1/A^ (here, A is the London 
penetration depth). From the temperature dependence 
of Hci, important information can be gained, partic- 
ularly, regarding the symmetry of the superconducting 
state, since the appearance of gap nodes strongly modifies 
the T-dependence of the superfluid density and thereby 
the penetration depth A(T). 

The experimental determination of Hd has been a 
challenging and controversial problem since the begin- 
ning of HTSC research. Not only the values reported 
for the first fiux penetration field Hp are scattered over 
a wide field range, but in strongly layered superconduc- 
tors with H perpendicular to the planes, striking fea- 
tures have been observed, such as a marked upturn of Hp 
for temperatures T < Tc/2 |]-||. The origin of the 
positive curvature of Hp at low temperatures has been 
the matter of various speculations. Different mecha- 
nisms have been proposed: Bean-Livingston surface bar- 
riers bulk pinning |||,^, a modification of the 
character of the field penetration in layered structures , 
a low-temperature enhancement of the superconducting 
order parameter in the normal layers 0, etc. 

For type II superconductors, there are at least two 
kinds of barriers which hinder the system from reach- 
ing a thermodynamic equilibrium state: (i) surface and 
geometrical barriers [p|-p^ governing vortex penetration 
into the superconductor, (ii) bulk pinning barriers gov- 
erning vortex motion in the superconductor. For the de- 
termination of Hci, the relevant barriers are those which 
govern vortex entry into the sample. Of particular in- 
terest in this context is the phenomenon of vortex creep 
over surface and geometrical barriers. From the theoret- 



ical point of view this subject has been discussed in Refs. 
pOf and however, to our knowledge no systematic 
experimental study has been carried out so far. 

In this letter we investigate the dependence of Hp 
on the magnetic field sweep rate dH/dt of isother- 
mal magnetization curves for the strongly layered 
Bi2Sr2CaCu208+5 (Bi2212) material. This has been 
done for crystals with different cross-sections along the 
c-axis. For specimens with ellipsoidal cross-sections geo- 
metrical barriers can be neglected [ p^|p^ , so that the rel- 
evant barriers for vortex entry are expected to be surface 
barriers. Indeed, for our specimen with an ellipsoidal- 
like cross-section, the Hp vs. dH/dt dependence ob- 
tained at high sweep rates is well described by vortex 
creep over surface barriers However, for decreas- 

ing dH/dt rates, these barriers are observed to collapse. 
The subsequent saturation of Hp at lower rates indicates 
that the system is in equilibrium with respect to surface 
barriers. This is interpreted as strong evidence that we 
have reached the lower critical field Hd in our measure- 
ments. The obtained low-temperature dependence of Hd 
is in good agreement with recent microwave absorption 
measurements of the penetration depth A On 
the other hand, for specimens with rectangular cross- 
sections, no significant vortex creep over the relevant bar- 
riers for vortex entry could be observed, in consistency 
with theoretical results for geometrical barriers . The 
low-temperature upturns of Hp observed for specimens of 
either cross-section are then explained in terms of mea- 
surements where the system is out of equilibrium with 
respect to barriers for vortex entry. 

The crystals under investigation have been grown with 
different techniques and have different shapes. The spec- 
imen with an ellipsoidal-like cross-section along the c- 
axis has been grown in a Zr02 crucible [T5| ] and is ap- 
proximately 1 X 1.3 X 0.05 mm'^ in size. The spec- 
imens with rectangular cross-sections along the c-axis 
have been grown with the traveling solvent floating zone 
method m& and are slightly larger in size. The experi- 
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ment is performed in a non-commercial SQUID magne- 
tometer where the sample is stationary in the pick-up 
coil. The field H is supplied by a superconducting coil 
working in a non-persistent mode. For the magnetiza- 
tion curve measurements, the sample is zero field cooled 
from above Tc and stabilized at a fixed temperature (the 
residual field of the cryostat is < lOmOe). Further, the 
field H is apphed at a fixed rate dH/dt. For Hp, we se- 
lect the field at which a deviation from Meissner shielding 
occurs (see inset of Fig. ph. 
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FIG. 1. For the specimen with an ellipsoidal-like 
cross-section, the first flux penetration field Hp vs. the applied 
magnetic field sweep rate dH/dt is displayed for different tem- 
peratures. The dotted lines are guides to the eyes. The data 
are scaled with the demagnetization factor A'^ = 0.96(±15%). 

Fig. shows the dependence of Hp on the field sweep 
rate dH/dt at three characteristic temperatures for the 
specimen with an ellipsoidal-like cross-section. A sharp 
step in the Hp vs. dH/dt curves is observed at all three 
temperatures. At high sweep rates, the curves display 
a finite slope which is most pronounced at T — 11 K 
and 17 K and to a lesser degree at T — 61 K. These 
slopes are in good agreement with the theoretical results 
of Burlachkov et al. [ pl| for vortex creep over surface bar- 
riers (see the analysis below). Proceeding from the creep 
regime towards lower rates, the sharp drop signals a col- 
lapse of the surface barriers (or equivalently a heat pulse) 
which leads to a rapid fiux entry, possibly in terms of flux 
jumps. Indeed, an avalanche type of fiux penetration into 
a type II superconductor has been reported by Duran et 
al. [0. On decreasing the sweep rates even further, the 
fiux penetrates into the sample through the occurrence of 
rare events when the applied magnetic field reaches Hd- 
We thus identify the saturated value of Hp with the lower 
critical field Hd and use it to obtain information about 
the nature of the superconducting state. 

Analogous measurements have been done on thin spec- 
imens with rectangular cross-sections. According to Zel- 
dov et al. |1U, for such specimens, geometrical barriers 



build up over a distance s from the sample edges, where s 
is the sample thickness. The energy required to overcome 
such a barrier is macroscopic ~ egs so that vortex creep 
over geometrical barriers is expected to be very weak 
(here Eq = (<i>o/47rA)^ and $o is the unit flux). As a mat- 
ter of fact, we were not able to detect significant vortex 
creep over geometrical barriers within the experimental 
time scale 10"^ < dH/dt < lO'^ Oe/s. 

We proceed with a brief discussion of vortex creep over 
surface barriers. As shown by Burlachkov et al. [ pT| , 
the surface barrier for pancake-vortices is given by t/ ~ 
£od\n {0.76Hc/ H), where He is the thermodynamic criti- 
cal field and d the interlayer distance. During the time t, 
thermal creep allows to overcome barriers of size U{t) c± 
T\n{t/to), where is a "microscopic" time scale [ [l8[ . 
Equating the two expressions one obtains the time de- 
pendence of the penetration field for the pancake-vortex 
regime, 

Hp{t) ~ H, {t/tor^^'"^ . (1) 

With the definition To = end/ \n{t/tQ) Eq. (|l|) becomes 
Hp ~ exp (— T/Tq). At higher temperatures creep 
proceeds via half loop excitations of vortex-lines [pT| . 
The creation of a half loop saddle configuration involves 
an energy U ~ TreegC In^ (jo/j)/2$oJ, where jo is the 
depairing current density and e = {m/M)^/'^ < 1 the 
anisotropy parameter. With a similar analysis as above, 
the time dependence of Hp for the vortex-line regime 
takes the form 

££qS, , ^ f Hr 



Hp{t) ~ He 



In" 



(2) 



2^/2 Tint /to 

where ^ is the coherence length. According to Ref . , 
this expression gives a temperature dependence Hp cx 
(Tc - Tf/^/T. The crossover between the two regimes 
occurs when the size of the half loop excitation is of the 
order of the interlayer distance d. This is expected to 
occur at a temperature T* w To(r*) In {d/eS). 

According to (Q), creep of pancake- vortices over sur- 
face barriers results in a linear dependence of \n{Hc/ Hp) 
on \n{t/to) with a slope T/eod. This dependence is given 
in Fig.g(a) using the data at high cycling rates in Fig.^. 
For T = IIK and 17K, the values of T/eod obtained 
from the fits are in good agreement with the values es- 
timated from the penetration depth Xab{T) (see Table |). 



TABLE I. T/eod values calculated with the penetration 
depth \ab{T) as obtained with formula Hd — ($o/4'7rA^j,)-ln k 
from the Hd data in Fig. ^ and values of T/eod obtained from 
the fits in Fig.&a). 





T= IIK 


T= 17K 


T/eod (calculated) 


1/34 


1/21 


T/eod (from fit) 


1/43 ±30% 


1/22 ±30% 
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On the other hand, for T = 61 K the above analysis for 
pancake-vortices is not satisfactory since the T/sod value 
obtained from the fit is one order of magnitude smaller 
than the calculated value. According to (||), for creep of 
vortex-lines over surface barriers {He/ Hp) ■ \n^{Hc/ Hp) 
vs. \n{t/to) is linear with a slope 2^/2T/neeo(,. For the 
T = 61 K data at high cycling rates in Fig. |l], such a rep- 
resentation is given in Fig.||(b). The slope 2-\/2T/7reeo^ 
obtained from the fit is 15 ± 3, in good agreement with 
the value 16, estimated with the help of the GL formula 
for the coherence length £,{T), with ^(0) ~ 25 A and 
an anisotropy parameter e ~ 1/100. From these con- 
siderations we conclude that the behavior of Hp at high 
sweep rates is determined by creep of pancake- vortices (at 
low temperatures) and vortex-halfloops (at higher tem- 
peratures) over surface barriers. The estimated value of 
the activation barrier for vortex entry is J7 ~ 300 K for 
temperatures T between 10 and 20 K and U ~ 1200 K 
for T ~ 60 K (here we set the time origin t = when 
the applied field H reaches Hci and assume to ~ 10^^ s). 
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FIG. 2. (a) T = IIK and 17K data of the upper plateau 
in Fig.|l] in a different representation. Here we choose 
the time origin t — nt H = Hd- The critical field 
Hc{T) ~ Hc{0) (1 - T/Tc) is calculated with Aab(O) as ob- 
tained from the Hd data in Fig.^ and the coherence length 
C(0) ~ 25 A. We assumed to ~ 10"^ s. (b) T = 61 K data at 
high sweep rates in Fig. ^ in a convenient representation. The 
lines in Figs. |(a) and (b) are fits according to Eqs. (jl|) and 
(^), respectively. 

In Fig.|^ we make use of the curves in Fig. |l] to deter- 
mine Hp{T) with three different criteria: (i) Hp data (o) 
are taken at low sweep rates so as to guarantee that they 
lie in the regime where the system is in equilibrium with 
respect to surface barriers (these data represent Hci{T)), 
(ii) Hp data (*) are taken at high rates so that they al- 
ways lie in the creep regime, and (iii) the field H is swept 
at the constant rate dH/dt = 8 • 10~'^Oe/s for all tem- 
peratures. In this case, the values of Hp (v) lie on the 



ffci-curve for T > 60 K, they go through the step-like 
crossover for temperatures 30 K < T < 60 K, and finally, 
for T < 30K they lie in the creep regime. From a com- 
parison of the curves in Fig.|^, it follows that the low- 
temperature upturn of Hp results from measurements 
where the system is out of equilibrium with respect to 
surface barriers. 
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FIG. 3. Temperature dependence of the first fiux 
penetration field Hp of the specimen with an ellip- 
soidal-like cross-section for different field sweep rates: 
(o) dH/dt < 1 ■ 10" 
(V) dH/dt = 8-10 
from the Meissner slope for T = IIK and 61 K. The upper 
curves (■) are measured at a rate dH/dt > 5 • 10~^ Oe/s and 
the lower ones (□) at dH/dt < 1 ■ 10"* Oe/s. The data are 
scaled with the demagnetization factor A*" = 0.96(±15%). 

As shown in Fig. ^ for increasing temperatures the T- 
dependence of Hp in the creep regime, (*) and (v) for 
T < 30K and (*) for T > 30K, undergoes a crossover 
from an exponential behavior to a weak power-law. This 
is in agreement with the results presented above and in- 
dicates that the crossover from creep of pancake- vortices 
to creep of vortex-halfloops over surface barriers takes 
place at T* w 30 — 40 K (in consistency with the esti- 
mate T* w To(T*)ln(d/e^) « 40K obtained using the 
parameters for Bi2212 and Int/to w 20). 

Making use of the correct Hd data for H _L ab, (o) 
in Fig.|^, we have investigated the temperature depen- 
dence of the penetration depth Xab using the expression 
Hci = ('i'o/47rA^j,) • luK, (here we assume that Iuk is 
T-independent). For IIK < T < 30 K, we find a lin- 
ear behavior XabiT) — Aah(O) cx T with slope dX/dT ~ 
lOA/K, where a linear extrapolation to the data gives 
a T = penetration depth XabiO) ^ 2700 A. Our 
measurements thus confirm the linear low-temperature 
dependence of Xab well as the slope dX/dT, which 
have been recently determined with microwave absorp- 
tion techniques on clean Bi2212 crystals ||l^,Q. A linear 
low-temperature dependence of Xab, as first observed by 
Hardy et al. [p9[ in YBCO crystals, is consistent with a 
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FIG. 4. Temperature dependence of Hp {H \\ ab) for a spec- 
imen with rectangular cross-section and size 2 x 3.9 x 0.05 mm^ 
for different sweep rates: (*) dH/dt > 6 • 10~^ Oe/s and (o) 
dH/dt < 6 • 10"^ Oe/s. The inset shows the Hp vs. dH/dt 
dependence at T = 12 K and 81 K. 

superconducting state with line nodes on the Fermi sur- 
face pO|, e.g., a d-wave symmetry of the order parameter. 

Similar investigations of Hp have been done for H \\ ab. 
The misalignment angle between H and the a6-planes 
has been estimated to be smaller than 2° from measure- 
ments at magnetic fields in the range -f/py <^ H < Hp±, 
(see Ref. here iJ^n and Hp± are the parallel and 

perpendicular penetration fields). For this configuration 
geometrical barriers are irrelevant. As shown in Fig.^, 
for temperatures T between 10 and 70 K, Hci{T) has 
an approximately linear behavior. According to Ref. [p2[ , 
for temperatures not so close to Tc, the correct de- 
scription of Hci in strongly layered materials is Hci = 
$o/(47rAa6Ac) • [\n{Xab/d,) + 1.12]. With this formula we 
calculated the penetration depth in c-direction Ac using 
the previously determined Xab data. For T < 40 K, Ac 
is approximately linear in T with a slope of 0.1/xm/K. 
The extrapolated T = value is Ac(0) ~ 15 ^m. Fig.| 
further shows Hp data at high sweep rates (*). Contrary 
to the case for H _L ab, no upturn of Hp is observed here 
for temperatures T < Tc/2. We attribute this difference 
to the absence of the strong pancake vortex creep regime 
at low temperatures for H \\ ab. From the Hp vs. dH/dt 
dependence in the inset of Fig.^, we obtain an activa- 
tion barrier for vortex entry U ~ 200 K at T = 12K and 
U ~ 1300 K at T = 81 K. 

Close to Tc a downward bending of iJci is observed 
for H ± ab a.s well as H || ab (see Figs.|| and This 
can be understood on the basis of an entropic downward 
renormalization of the vortex line free energy due to fluc- 
tuations of the order parameter around its mean-field 
form . The decrease in the free energy fi — £i — Tsi 
then leads to a drop in Hd as T approaches Tc (here, £; 
is the line energy of the vortex excitation and si is the 



line entropy). 

To conclude, isothermal magnetization curves were 
measured on Bi2212 crystals for configurations of neg- 
ligible geometrical barriers. At very low sweep rates we 
found saturated values of the penetration field Hp which 
we interpret as the true lower critical field Hd . The low- 
temperature dependence of Hd is consistent with recent 
magnetic penetration depth measurements |l^,|ljl sug- 
gesting a superconducting state with nodes in the gap 
function. Based on the results obtained for Bi2212, we 
argue that the frequently reported low-temperature up- 
turns of Hp in strongly layered superconductors with H 
perpendicular to the planes [0-|6| can be explained in 
terms of measurements where the system is out of equi- 
librium with respect to barriers for vortex penetration. 
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